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Thermodynamic characteristics of Fermi systems are investigated in the vicinity of a phase tran- 
sition where the effective mass diverges and the single-particle spectrum becomes flat. It is demon- 
strated that at extremely low temperatures T, the flattening of the spectrum is reflected in non- 
Fermi-liquid behavior of the inverse susceptibility x" 1 ^) ~ T a and the specific heat C{T)/T ~ 
T~ Q , with the critical index a = 2/3. In the presence of an external static magnetic field H, both 
these quantities are found to exhibit a scaling behavior, e.g. x~ X {T,H) = x~ X (T, 0) +T 2/3 F(H/T), 
in agreement with available experimental data. 
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Manifestations of non- Fermi- liquid (NFL) behavior ob- 
served at extremely low temperatures T in a number 
of strongly correlated Fermi systemsr'S' 3 *^' 6 *^'^' 9 ^ no- 
tably 3 He films and heavy-fermion metals, provide valu- 
able clues to a fundamental microscopic understanding 
of these systems. This anomalous behavior is commonly 
attributed to antiferromagnetic spin fluctuations, and nu- 
merous strategies have been advanced (e.g. Ref.lll|) based 
on integrating out all degrees of freedom except spin fluc- 
tuations. However, the spin- fluctuation model (SFM) is 
unable to reproduce the results of precise measurements 
of the inverse spin susceptibility, which, on the "metallic" 
side of the critical point, is found to obey v~ 1 (T) ~ T a 
with a critical index a < 1 (cf. Refs. 19 12). Difficul- 
ties are likewise encountered for the thermal expansion 
dataii The SFM also fails to explain 6 : 9 : 10 : 13 the scaling 
exhibited by x(T, H) in external, sometimes tiny, static 
magnetic fields H . Moreover, the SFM falls short when 
confronted with other generic features revealed by exper- 
iment, most notably a divergence of the effective mass 
M* (p) at a critical density p = p^ and the emergence of 
NFL features even before the critical poin L 1 : 3 : 4 : 5 

We are therefore compelled to explore a different strat- 
egy, in which NFL anomalies are attributed to fermionic 
degrees of freedom and specifically associated with flat- 
tening of the single-particle spectrum £(p) in the imme- 
diate vicinity of Poo- Near this point, the FL formula 
£fl(p) = Pf{p — Pf)/M*(p) must be supplemented by 
terms nonlinear in p — pp, since £fl(p; p) vanishes iden- 
tically at poo. Here we shall study consequences of the 
alteration of on the "metallic" side of the phase tran- 
sition. Importantly, as argued below, the ratio of the 
damping 7(e) of single-particle excitations to relevant en- 
ergies e remains relatively small in this regime. Accord- 
ingly, the Landau quasiparticle formalism still applies. 

We focus our attention on the real part of the AC 
spin susceptibility x(T,uj — > 0) given by the familiar FL 
formulaic 



X(T, p) = Xo(T, p)/[\- 5 o(p)n (T, p)] 



(1) 



where go is the zeroth harmonic of the Landau spin-spin 
interaction xo = — P%^o, & n d H-b is the Bohr magne- 
ton. To be definite, we treat the homogeneous three- 
dimensional (3D) case, with 



n (T) 



dn{(j(p)) = pp 



n(0[l-n(£)]|k: 



(2) 

where n(£(p)) = 1/[1 + exp(£(p)/T)] 1 is the quasiparti- 
cle momentum distribution and dv = 2d 3 p/ (2ir) 3 . While 
Eqs. and J2J) are formally identical to familiar text- 
book formulas, the function Hq(T) can differ profoundly 
from its FL realization due to complicated character of 
the group velocity d£(p)/dp in the momentum region 
where \t;{p)\ — T, which is dominant in the integral (J2J. 
To proceed, we invoke the well-known expression^ 
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M 



/(P,Pi) 5 dv i> ( 3 ) 
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connecting the spectrum £(p) and the momentum dis- 
tribution n(£) through the Landau interaction function 
/(p,Pi). At T = 0, Eq. © implies that the effective 
mass M* is related to the first harmonic /i(pf,Pf) of 
the interaction function by 



M/M*{p,T = 0) = 1 - F°(p)/3 = D(p) . 



(4) 



in which Ff{p) = fi(pF,PF)N with N = ppM/n 2 . 
Thus M*(p,T — 0) diverges at the critical density 
Poo, where D(poo) = 0, while at nonzero temperatures, 
M*(poo,T) already has a finite valued Its T dependence 
is found by expanding relevant quantities on both sides 
of Eq. J2J in Taylor series, thereby obtaining 



d^/dp ~ p F /M*(p,T) + v 2 { P -pf) 2 /M Pi 



(5) 



with 
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(f[p 2 F + 2f lPF )(p- PF ) 



2 
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Here v 2 = -p 3 F Mf{' /6tt 2 , f[ = [df 1 (p,p F )/dp} p=pF , and 
fl' = [d 2 h{p,PF)/dp 2 ] p=pF . The term Vi(p)(p - p F ), 
which provides the contribution v\{p)(p — p F ) 2 /2 to 
has been dropped in writing Eq. J5J), since i>i(poo) must 
vanish; otherwise the function £(p, poo) has the same sign 
below and above the Fermi surface, and the Landau state 
becomes unstable before p reaches poo- 

Eq. Q can be simplified using particle- number conser- 
vation, expressed approximately as 



/ \Pf(P Pf) + (P- Pf) 2 } ^^dp = . 
Inserting this relation into Eq. JSJ, we have 

M nf , f ( P -PF) 2 dn(Z(p)) , 
D(p)-v 2 / — — dp. 



M*(T,p) 



Pf 



dp 



(7) 
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In terms of the new variables r = 2TM/pp, x = £(p)/T 
and y = (p - p F )(w 2 /3MTp F ) 1 / 3 , Eq. © reduces to 
M/M*(T,p) = D(p) + var 2 / 3 , while Eq.© becomes 



d£{p, T, p) _ p£ 
dp M 



D{p) + VT 2 / 3 [a + y 2 )] 



which is recast to a relation between x and y, 



x = y 



Zv- 1 D( P )t- 2/3 + 3a + y 2 



where v = (9V2/4) 1 / 3 and 



y 2 (x)e x (l + e x )- 2 dx. 



(9) 



(10) 



(11) 



Having solved Eqs. I|10[l and (|11|) . the polarization oper- 
ator Ho = — NqPq is straightforwardly evaluated to yield 

X (T,p)=fx 2 B N P (T,p)/[l + G (p)D(p)P (T,p)} , (12) 

where G = g p F M*(T = 0)/tt 2 , and 

e x (l + e x )~ 2 



Po(T,p) 



[a(T,p)+y 2 (x)]vT 2 /3 +D(p) 



dx . (13) 



In what follows we address systems without ferromag- 
netic ordering where 1 + Go > 0. In this case, Po(T,p) 
crucially depends on the ratio T/\p — Pool- When T 
drops to while \p — poc \ is fixed, one obtains Pq(T, p) — 
D-\p), and hence X (T = 0,p ± Pao ) = p 2 B N /[D(p)(l- 



Go)] ~ \ P 



In the opposite case, D{p oc ) = 



and calculations yield M/M*(T, p^) ~ a(poo)^r 2 / 3 , and 
Po(T 7Poo ) = bip^v- 1 ^ 2 ' 3 with a{ Poo ) = 0.5 and 
b{Poo) — 1-2, so that 

X (T, Poo ) = X o(T, poo) = UplN^T- 2 ' 3 . (14) 



This implies that the critical index a specifying the low-T 
dependence of the inverse spin susceptibility x -1 (T, poo) 
is 2/3. 

It can be verified that the corresponding ID and 2D 
equations are identical in form to Eqs. (|10f) and Ullfl . 
derived here for the 3D case. 

The above treatment can be extended to other ther- 
modynamic properties, e.g. the specific heat given by 



(15) 



Manipulations similar to those applied to Hq (T, p) yield 



C(T,p) 
T 



N 



x[x — 2ay(x))e x (l + e x ) 
[a + y 2 (x)]vT 2 / 3 +D{p) 



-dx . (16) 



Thus we see that at T — > the ratio C(T)/T possesses 
the same NFL behavior as the spin susceptibility. 

In electron systems of solids, as a rule there exist sev- 
eral bands that cross the Fermi surface, the effective 
mass diverging only in one or two of them. Other bands 
still contribute as Landau theory dictates providing ad- 
ditional T-independent terms in x(^) an d C(T)/T and 
modifying somewhat the prefactors in the singular parts. 
As a result, the behavior T -2 / 3 implied by Eq. ljTfll holds 
in dealing with the T-dependent excesses Ax{T, poo) and 
AC(T, poo)/T. One may expect the same to be true for 
anisotropic systems, e.g. for 2D liquid 3 He, where the 
enhancement of M* is associated with the crystal lattice 
of the substrate^ 

To ascertain the relevance of the basic result 
(|14fl to real strongly correlated Fermi systems, we 
make use of data&2iifi for the heavy-fermion met- 
als YbRh(Sio. 95660.05)2 and CeRu2Si2, where the T- 
independent FL terms are rather small. These data, as 
seen from the left panel of Fig. are well reproduced by 
the proposed model. 

In 2D liquid 3 He the FL term in x(^) is significant. 
Accordingly, in the right panel of Fig. ^ we plot results 
for the excess Ax{T) at the two densities p = 0.036^4~ 2 
and p = 0.052A~ 2 . In doing so, we suggest that the flat 
portion in the spectrum £(p) forms in certain regions 
adjacent to the Fermi line of 3 He, and these regions pro- 
duce contributions A%(T) ~ T~ 2 / 3 , while the remainder 
generates only the FL terms. We suggest that with in- 
creasing density the latter region shrinks and vanishes at 
the point where liquid 3 He solidifies. Our calculations ar 
compared with experimental data— in the right panel of 
Fig.IU 

Let us now examine the role of damping effects. When 
evaluating thermodynamic properties, characteristic en- 
ergies e are of order of T. We estimate the damping 
7(5 ~ T) with the help of the standard FLT formula 



7 ( £ ~T)~|r| : 



[M*(T)] 3 T 2 



(17) 



containing the interaction amplitude T and the effective 
mass M*, which according to Eq. @, depends on T as 



3 





FIG. 1: Left panel, bottom-left axes: Inverse magnetic suscep- 
tibility of YbRh(Sio.95Geo.Q5)2 as a function of temperature. 
Experimental data of Ref. [H are shown as short dashes; the 
solid curve is the current prediction. Top-right axes: tem- 
perature dependence of magnetic susceptibility of CeRu2Si2 
in a magnetic field of 0.02 mT. Experimental data of Ref. ^3 
are denoted by circles, and the theoretical prediction by the 
dashed curve. Right panel: Spin susceptibility excess, divided 
by the FL contribution to the spin susceptibility, evaluated for 
3 He films at two densities (indicated near the curves). Exper- 
imental data from Ref. |2j appear as solid and open circles, 
while solid curves trace the predictions of the current theory 
at low T. 



T~ 2 / 3 . It is important to recognize that the customary 
replacement of T by the bare interaction V, although 
legitimate in ordinary Fermi liquids, is erroneous in the 
limit of strong correlations. The source of the error is 
the huge enhancement of the density of states, which sup- 
presses r and makes its magnitude quite insensitive to the 
bare interaction. 18 To wit: summation of ladder diagrams 
in the particle-hole channel gives T = V/[l + VN(0)], 
where N(0) ~ ppM* is the density of states. In the limit 
N(0)V > 1, one obtains T ~ iV _1 (0). Inserting this 
result into Eq. I|17|l . we obtain 



y(e ~T)/T ~TM*(T)/p 



-1/3 



< 1 



(18) 



affirming the applicability of FLT theory to our problem. 

Imposition of a static external magnetic field H 
brings into play a new dimensionless parameter R — 
HbH/T and opens another arena for testing the model. 
The function n(^(p)) entering Eq. is then re- 

placed by [n(£ + (p)) + n(£_0)]/2, where n{£_±(p)) = 
[1 + cxp(£(p)/T ± R/2)]~ l . In turn, is determined 
from Eq. J3J), the effective mass being the same for both 
spin directions at sufficiently weak H. Proceeding as be- 
fore, we find 



M/M*(T,H, Poo ) = VT 2 ^a{R), 



where 
a{R) 



y 2 (x) 



^x+R/2 



,x-R/2 



(19) 



dx . 



(l + e x + R / 2 ) 2 (l + e x - R / 2 ) 2 

In the limit T — > or equivalently R — > oo, the so- 
lution of this equation takes the analytic form M*(T — 
0, H, p l30 ) ~ H 2 / 3 , confirming a result of Ref.lT^ Thus at 
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FIG. 2: Top panels: Normalized magnetic susceptibility 
x(T,H)/x(Tp) (top-left panel) and normalized magnetiza- 
tion M(T, H)/M(Tp) (top-right panel) for CeRu2Si2 in mag- 
netic fields 0.20 mT (squares), 0.39 mT (triangles), and 0.94 
mT (circles), plotted against normalized temperature T/Tp 
(Ref. IT??) , where Tp is the temperature at peak susceptibil- 
ity. The solid curves trace the universal behavior predicted 
by the present theory. Bottom panel: The normalized ra- 
tio C(T, H)T M /C(T M )T for YbRh(Sio. 95 Geo.o5)2 in magnetic 
fields 0.05 T (squares), 0.1 T (triangles), 0.2 T (circles) versus 
the normalized temperature T/Tm (Ref.Qlth where Tm is the 
temperature at maximum ratio C(T, H) /T . The solid curve 
shows the prediction of our theory. 



sufficiently low temperatures, imposition of static mag- 
netic fields satisfying /j,bH > T renders the effective mass 
M*(T, H, poo) finite, promoting the recovery of the Lan- 
dau theory. Along the same lines, at the critical density 
Poo we may establish that the magnetic moment and AC 
spin susceptibility display a scaling behavior, e.g. 

X ac(T, H, poo) = M |iVo^- 1 6(i?)T- 2 / 3 , (21) 



where 



x+R/2 



-R/2 



(\J re x+R/2\2 (i +e x-R/2y 



dx 



a{R)+y 2 {x) 

Turning to the specific heat, and following the same 
path as taken above for x{T, H) at the critical point, we 
are led to the expression 

C(T, H, poo)/T = No^ciR)^ 3 , 



(23) 



where 



c(R) 



[{x+R/2) 2 -2{x+R/2)a{R)y{x)}e 



x+R/2 



(1 + e x + R / 2 ) 2 
{x-R/2) 2 -(x-R/2)a{R)y{x)]e x - R l T 

( 1 + e x-fl/2)2 



dx 



a(R)+y 2 (x) " 
(24) 
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A scaling behavior of this kind has been discussed by 
Coleman and collaborators*^^ 

For finite H, the curve describing xac(T,H) acquires 
a maximum at some temperature Tp, as does the curve 
for C(T,H)/T, at some temperature Tm- This behavior 
is associated with the suppression of the divergent NFL 
terms - T~ 2 / 3 in x{T,H) and C{T,H)/T and recov- 
ery of the FL behavior at static magnetic fields in which 
the Zeeman energy splitting psH exceeds T. Follow- 
ing Ref. we present in Fig. 2 the results of numerical 
calculations of these quantities as functions of the nor- 
malized temperatures T/Tp and T /Tm to demonstrate 
that our model reproduces the experimental scaling be- 
havior of the spin susceptibility^ of the heavy fermion 
metal CeRu2Si2 and that of the specific heals^ of the 
heavy fermion compound YbRh(Sio. 95660.05)2, without 
adjustable parameters. It should be emphasized that the 
curves given in Fig. 2 remain the same whether one is 
working with heavy fermion metals or 2D liquid 3 Hc. 
This universality predicted by our model, can be tested 
with the aid of a new apparatus^ designed for measure- 
ments of thermodynamic properties of 2D liquid 3 He in 
static magnetic fields. 

It is worth noting that, as seen from Eq.J2J), on the 
metallic side of the phase transition the NFL behavior 
(|14fl holds in a density domain Ap ~ Poor 2 / 3, adjacent 
to the critical point. On the other hand, as we plan to 
demonstrate in our next article, on the insulating side 
of the phase transition, the range of the domain where 
the NFL term ~ T~ 2 / 3 prevails is substantially larger, 
Ap ~ poor 1 / 3 . 

Let us compare our results with those available within 
the antiferromagnetic scenario First, the proposed flat- 
tening mechanism for NFL behavior adequately explains 
the low-T data on the spin susceptibility producing 
X _1 (T) ~ T a with a ~ 2/3 in the critical density re- 
gion, while the SFM fails to provide a < 1. Second, 
the flattening mechanism explains the scaling behavior 
X -!(T) ~ T a F(H/T) of the spin susceptibility in static 
magnetic fields, whereas such a scaling property does not 
arise in the SFM. Third, within the model developed 
here, FL behavior can be recovered at low T close the 
critical point by imposing a tiny magnetic field satisfy- 



ing pbH > T. In the SFM there is no such provision for 
reinstating Fermi-liquid theory. 

Finally, we discuss the relevance of ferromagnetic fluc- 
tuations to the NFL behavior, claimed in Rcfs. IglfTol Ex- 
perimental data for most heavy-fermion systems show 
no evidence of ferromagnetism. The same is true for 
2D liquid 3 Hei 3 « This implies that close to the critical 
point, the respective Pomeranchuck stability condition 15 
1 + Gq(p) > is not violated, in spite of the diver- 
gence of the effective mass at p = p^. The value 
of Go can be estimated from the Landau sum rulei^ 
J2 L {F L /[1+F L /(2L+1)}+G L /[1+G L /(2L+1)}} = by 
assuming all the harmonics Fl to be finite, apart from 
F X =F^M* /M~3M* /M. The Landau sum rule is then 
satisfied provided Go— — 0.75. This value of Go gives rise 
to a marked enhancement of the Sommerfeld ratio in the 
case pbH > T, to a around 12 that changes little if the 
zcroth harmonic Fq is taken into account. Even so, such 
an enhancement does not in itself ensure the relevance 
of ferromagnetic fluctuations to the NFL behavior. In- 
deed, dense 3D liquid 3 He, in which Go — —0.8, shows 
no deviations from FLT predictions as T — > 0. 

In summary, we have studied the impact of flattening 
of the single-particle spectrum on the spin susceptibility 
x(T) and the specific heat G(T) of strongly correlated 
Fermi systems. When the density approaches the critical 
value at which effective mass becomes infinite, Fermi- 
liquid theory progressively fails as the deviant compo- 
nents in x(T) and G(T) grow to dominance. We have 
explicated the implied critical behavior of these compo- 
nents and derived a universal relation between them. Fi- 
nally, our analysis has revealed the scaling behavior of 
x(T,H) and C(T, H)/T in the presence of an external 
magnetic field H. Numerical results based on this theo- 
retical picture are in agreement with experimental data 
on 2D liquid 3 He and several heavy-fermion metals. 
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